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In this paper we consider a superconducting film modeled 
by the Ginzburg-Landau model, confined between two paral- 
lel planes a distance L apart from one another. Our approach 
is based on the Gaussian effective potential in the transverse 
unitarity gauge, which allows to treat gauge contributions in a 
compact form. Using techniques from dimensional and zeta- 
function regularizations, modified by the confinement condi- 
tions, we investigate the critical temperature as a function of 
the film thickness L. The contributions from the scalar self- 
interaction and from the gauge fluctuations are clearly iden- 
tified. The model suggests the existence of a minimal critical 
thickness below which superconductivity is suppressed. 

PACS number(s): 74.20.-z, 05.I0Cc, 11.25.Hf 



I. INTRODUCTION 

It is currently assumed to be a good approximation to 
neglect magnetic thermal fluctuations in the Ginzburg- 
Landau (GL) model, to explore general features of su- 
perconducting transitions. However, this approximation 
excludes the study of the so called charged phase transi- 
tions. They only can be investigated when fluctuations 
of the gauge field are taken into account, which make 
explicitly appear in the thermodynamic quantities the 
coupling constant for the interaction between the scalar 
field and the gauge field (the charge of the boson) . This 
is a hard problem to be directly faced and many attempts 
have been done to go beyond the pioneering work of ref. 
[1]. In what concerns physical situations already present 
in the literature, a large amount of work has been done 
on the Ginzburg-Landau model applied to the study of 
superconductors, both in the single component and in the 
N -component versions of the model, using the renormal- 
ization group approach. The interested reader can find 
an account on the state of the subject for both type-I 
and type-II superconductors and related topics in Refs. 
[2-8]. 

In this paper, we intend to take into account gauge fluc- 
tuations in a study of superconducting films. We consider 
the Ginzburg-Landau model, the system being submit- 
ted to the constraint of confinement between two parallel 
planes a distance L apart from one another. From a phys- 
ical point of view, for dimension d = 3 and introducing 
temperature by means of the mass term in the Hamil- 
tonian, this corresponds to a film-like material. We in- 
vestigate the behaviour of the system taking into account 
gauge fluctuations, which means that charged transitions 
are included in our work. We are particularly interested 



in the problem of how the critical behaviour depends on 
the film thickness L. This study is done by means of the 
Gaussian Effective Potential (GEP) as developed in Refs. 
[9-12], together with a spatial compactification mecha- 
nism introduced in recent publications [13,14]. 

Effects of spatial boundaries on the behaviour of phys- 
ical systems appears in several forms in the literature. 
For instance, there are systems that present defects, as 
domain walls created for instance in the process of crys- 
tal growth by some prepared circumstances. At the level 
of effective field theories, in many cases this can be mod- 
eled by considering a Dirac fermionic field whose mass 
changes sign as it crosses the defect, what means that 
the domain wall can be interpreted as a kind of a criti- 
cal boundary [15,16]. Questions concerning stability and 
the existence of phase transitions may also be raised if we 
enquire about the behaviour of field theories as function 
of spatial boundaries. The existence of phase transitions 
would be in this case also associated to some spatial pa- 
rameters describing the breaking of translational invari- 
ance, in our case the distance L between planes confining 
the system (a superconducting film of thickness L). In 
particular the question of how the superconducting criti- 
cal temperature could depend on the thickness of the film 
can be raised. 

In the next Section, we apply the functional ap- 
proach of the Gaussian effective potential formalism to 
the Ginzburg-Landau model, obtaining the mass, which 
obeys a generalized "Gaussian" Dyson-Schwinger equa- 
tion. In Section III, we extend the formalism of Sec- 
tion II to the GL model, confined between two parallel 
planes and a study of the critical behaviour of the sys- 
tem is performed. In particular, we obtain an expression 
for the critical temperature as a function of of the spac- 
ing between the planes (the film thickness). Finally, we 
summarize our results in Section IV and a qualitative 
comparison with some experimental observations is done. 



II. THE GAUSSIAN EFFECTIVE POTENTIAL 
FOR THE GINZBURG-LANDAU MODEL 

We begin by briefly presenting the study of the U(l) 
Scalar Electrodynamics in the transverse unitarity gauge, 
along the lines developed in ref. [17]. We start from the 
Hamiltonian density of the GL model in Euclidean d- 
dimensional space written in the form [18], 

H' + \\{d ll -ieA^ 



+ -m^| 2 + A(|vI/| 2 ) 2 , 



(2.1) 



where 



where f is a complex field, and mo is the bare mass. 
The components of the transverse magnetic field, = 
d^Av — dvAp {fi,v — l,...,d) are related to the d- 
dimensional potential vector by the well known equation, 



^ V" = |V x A| 2 . 



(2.2) 



In order to obtain only physical degrees of freedom, we 
can introduce two real fields instead of the complex field 
'J, assuming a transverse unitarity gauge. We can define 
the field in terms of two real fields, as = 0e* 7 , to- 
gether with the gauge transformation A — ► A — 1/eW-y. 
The unitarity gauge makes the original transverse field 
to acquire a longitudinal component A^ proportional to 
V7. Then the original functional integration over \f and 
\E'* in the generating functional of correlation functions, 
becomes an integration over <p, A T and A L . The longitu- 
dinal component of the vector potential can be integrated 
out, leading to the generating functional (up to constant 
terms), 



D<j) DA T ]exp 



d d xH 



x j4> 



(2.3) 



where the Hamiltonian is 



+ ieW + i(VxA) 2 + l(V.A) 2 . 



(2.4) 



We have introduced above a gauge fixing term, the limit 
e — > being taken later on after the calculations have 
been done. In Eq.(2.4) and in what follows, unless ex- 
plicitly stated, A stands for the transverse gauge field. 

The Gaussian effective potential can be obtained from 
Eq.(2.4), performing a shift in the scalar field in the form 
(j) ~ (p + ip, which allows to write the Hamiltonian in the 
form 



H = Tto + Hi 



(2.5) 



Vi = m^ip + 4A93 3 , 

1 2 2 , r \ 2 1 r->2 

v 2 = ^"W + 6Xip - -il , 

v 3 = 4Xip, 
v 4 = X. 



(2.8) 
(2.9) 
(2.10) 

(2.11) 
(2.12) 



It is clear from Eqs. (2.5, (2.6) and (2.7), that H de- 
scribes two interacting fields, a real scalar field 4> of mass 
fl and a real vector gauge field A of mass A. 
The effective potential, which is defined by 



Veff[<p] 



1 

V 



-lnZ+ / d xjtp 



(2.13) 



where V is the total volume, can be obtained at first 
order from standard methods from perturbation theory. 
One can find, from Eqs. (2.3), (2.6) and (2.7), 



Veff[<fi] 



It (SI) + 2I d (A) + -m 2 ^ 2 + Xip 4 

+1 [to 2 + 12A<^ 2 + 6XI d (n)] i d (n) 
+ [e 2 ^ 2 +/ i (Q))-A 2 ]/ d (A), (2.14) 



where, 



d d k 



(2ir) d k 2 + M 2 ' 



(2.15) 



Id ^ = \j (^ ln (fc 2 + M 2 ), (2.16) 

with k = (fci, kd) being the (i-dimcnsional momentum. 

The Gaussian effective potential is derived by the re- 
quirement that V e ff[ip] must be stationary under varia- 
tions of the masses A and f2. This means that values 
and A for the masses Q and A should be found such that 
the conditions, 



with Ho being the free part and Hint the interaction part, 
given respectively by 



Hq = 
+ 

and 



i(W) 2 + i^ 2 



i(VxA) 2 + iA 2 ^ + l(V-A) 2 



4 1 

n=0 

+ ±e 2 4>A»A^+±e 2 A^<p 2 , 



(2.6) 



(2.7) 



dV f 



eff 



on 2 

dV eff 



0, 



dA 2 



(2.17) 
(2.18) 



be simultaneously satisfied. These conditions generate 
the gap equations, 



Q = ml 



UXip 2 + 12XI d (n) + 2e 2 / d (A), (2.19) 
~K = e 2 ^ 2 + e 2 I d (U). (2.20) 



Replacing and A in Eq.(2.14) by the solutions and 
A, of Eqs.(2.19) and (2.20) we obtain for the GEP the 
formal expression, 



Veff[<P] 



+ 2I d (A) + im^ 2 + \tp 4 
3X[I d (T})} 2 -e 2 I d (T})I d (A). 



m = O n - 



(2.21) 



Notice that Eqs.(2.19) and (2.20) are a pair of very in- 
volved coupled equations, and no analytical solution for 
them has been found, they can be solved only by numer- 
ical methods. Later on we will see that this difficulty, in 
the limit of criticality can be surmounted. 

Next we intend to write an expression for the Gaussian 
mass, m, obtained in our case from the standard prescrip- 
tion, as the second derivative of the Gaussian effective 
potential for tp = 0. To calculate the second derivative of 
V e ff with respect to tp, we remark from Eqs.(2.19) and 

2 2 

(2.20) that ft and A also depend on ip according to the 
relations 



(2.30) 



Hence, we see from the gap equation (2.27) that the 
Gaussian mass obeys a generalized "Gaussian" Dyson- 
Schwinger equation, 

m 2 = m 2 + I2\l d (m) + 2e 2 / d ^e 2 I d (m)^j . (2.31) 

This expression will be used later to describe the system 
in the neighbourhood of criticality. 



III. CRITICAL BEHAVIOUR OF THE 
CONFINED GINZBURG-LANDAU MODEL 

A. The effect of confinement 
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dA 2 
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1 + 6X1^(0.) 



2eV 



-e 2 I d 
2 " 



where 



it^M) = 2 



r d d k 
J W) 



r) d (k 2 + M 2 ) 2 ' 
Replacing Eq.(2.23) in (2.22) we get, 

dTf _ [24A-2e 4 /f 1 (A)] tp 

dtp ~ 1 + [6A - \e 4 I d 1 (tT)] 1^(11) ' 



(2.22) 
(2.23) 

(2.24) 
(2.25) 



and the second derivative of the GEP with respect to tp 
is given by, 



eff 



dtp 2 



= to 2 , + 12Xtp 2 

+ 12A/ d (O) + 2e 2 I d (K) + 2eV^-i(A) 
[6A + |e 4 Jf !(A)] [24A - 2e 4 /f 1 (A)] tp 2 



\+[s\-\e*i-x{n)]iii(n) 



(2.26) 



Thus we have the formula for the Gaussian mass, 

7 eff 



m 2 



d 2 V e 



dtp 2 



<p=0 



m 2 + l2\I d (n ) + 2e 2 I d {A ), 



(2.27) 



where On and Ao are respectively solutions for O and A 
at tp = 0, explicitly, 



Oq = m 2 + 12A/ d (O ) + 2e 2 / d (A ) 
A^-e 2 / d (n ). 



(2.28) 
(2.29) 



Therefore, from Eqs.(2.27) and (2.28) we get simply, 



Let us now consider the system confined between two 
parallel planes, normal to the x^-axis, a distance L 
apart from one another and use Cartesian coordinates 
r = (xd, z), where z is a (d— l)-dimensional vector, with 
corresponding momenta k = (fcd,q), q being a (d — 1) 
-dimensional vector in momenta space. In this case, the 
model is supposed to describe a superconducting mate- 
rial in the form of a film. Under these conditions the 
field (f>(xd,z) satisfies the condition of confinement along 
the x^-axis, tp(xd = 0, z) = <p(xd = L,z) = 0, and 
should have a mixed series-integral Fourier expansion of 
the form, 



oo „ 

)= Yl d d -V(q)e-™- iq - z £K,q) 



(Xd,Z 



(3.1) 

where uj n = 2im/L and the coefficients c n and 6(q) cor- 
respond respectively to the Fourier series representation 
over Xd and to the Fourier integral representation over 
the (d— l)-dimensional z-space. The above conditions of 
confinement of the x^-dependence of the field to a seg- 
ment of length L, allow us to proceed with respect to 
the Xd-coordinate, in a manner analogous as it is done 
in the imaginary-time Matsubara formalism in field the- 
ory. The Feynman rules should be modified following the 
prescription, 



dk d 
~2T 



1 

I 



2nn 

— = uj n . (3.2) 



We emphasize however, that here we are considering an 
Euclidean field theory in d purely spatial dimensions, 
we are not working in the framework of finite tempera- 
ture field theory. Temperature is introduced in the mass 
term of the Hamiltonian by means of the usual Ginzburg- 
Landau prescription. 

For our purposes we only need the calculation of the 
integral given in Eq.(2.15) in the situation of confine- 
ment of the present section. With the prescription (3.2), 



the equation corresponding to Eq.(2.15) for the confined 
system can be written in the form, 



1 



+ 00 

E 



4tt 2 L " j q' 

n— — 00 



q 



(3.3) 



where % = fci/27r, a = l/L 2 and c 2 = M 2 /47r 2 . 

Eq.(3.3) can be treated within the framework of the 
formalism developed in Refs. [13], [14]. Using a well 
known regularization formula [19], we can write Eq.(3.3) 
in the form 



3-d 



-,a 



(3.4) 



where Af (^2^, a) is one of the Epstein-Hurwitz zeta- 
functions, defined by [20] 



+00 



Ak{y;{ai}) 



J2 (om? + ... + a K n 2 K + c 2 ) v , 



ni,...,nK=—oo 



(3.5) 



with Re(v) > K/2 (in our case Re(d) < 2). The 
Epstein-Hurwitz zeia-function can be extended as a 
meromorphic function to the whole complex z/-plane (for 
us, to all values of the dimension d), and we obtain af- 
ter some rather long but straightforward manipulations 
described in detail in [13], the expression, 



I${M) = 2-* tt 1 



n=l 



K 



i (MLri) 



(3.6) 



where are the Bessel functions of third kind. 

B. Critical behaviour 

Let us take M = m in Eq.(3.6) and let us restrict our- 
selves to the neighbourhood of criticality, that is, to the 
region defined by m « 0. Then the asymptotic formula, 

K v {z)* l -T{v)(^y V , (z~0) (3.7) 

allows to write Eq.(3.6) in the form 

d\ 1 



/ o d (m«0) 



r 1-; 



2 1 L d ~ 2 



C(d-2), (3.8) 



where £(d — 2) is the Riemann zeta-function, £(d — 2) = 
Er=i(V" d " 2 ) 5 defined for d > 3. For d > 3, in the sense 
of the analytic continuation in dimension of the Espstcin- 
Hurwitz z eta- functions, we obtain the expression, 



7 d (m«0) 



1 1 



2^L 



C(d-2). 



(3-9) 



The integral ig (A = ^Je 2 lg{m)j , w hich enters 

Eq.(2.31), must be considered carefully. For a dimension 
d > 3, we get, 



7 d (A ) « 2-t7T-4 



n=l 



2^r 

Ao 
nL 



Ki(A Ln) 



(3.10) 



or, using the exact expression for the summation in the 
above equation, 

n=l K ' 

(3.11) 



Eq.(3.10) becomes (with A = ^e 2 lg(m)) 
I d a Ue 2 !^ ' 



20? 



-V2n- In (l 



(3.12) 



However, notice that if we are in the limit m w 0, we see 
replacing Iq(Wi) from Eq.(3.9) in the exponential, that 
the logarithm in the last term of Eq.(3.12) will disappear 
at d = 3, due to the divergence of ((d — 2) as d — > 3. 
Hence, Eq.(3.12) becomes simply, for a dimension d > 3, 



1 



2tt 1 /4^/2 is 



-CM^-2). (3.13) 



Thus we can write the Gaussian gap equation (2.31) in 
the neighbourhood of criticality in the form, 



m 2 



24 1 
m 2 + —\-C(d-2) 



1 



e a -^d(rf-2). 



(3.14) 



For m = 0, Eq.(3.14) defines a critical equation for d > 
3. But it is well known that the only singularity of the 
zeta-function Q(z) is a pole at z = 1, which would make 
Eq.(3.14) meaningless as it stands for d = 3, just the 
physically interesting situation. 

However, we can give a physical sense to Eq.(3.14) for 
d = 3, by means of a regularization procedure. This can 
be done using the formula, 



lim 



CO*) 



= 1, 



(3.15) 



where 7 is the Euler constant, to define for d > 3 a new, 
dressed mass m, related to the former bare mass by, 



2 2 
m = m — 



24A 



V^i(d - 3) 



+ 



p=i 



YcW^ 



(d - 3)p ' 



(3.16) 



where C? are appropriate generalizations of the coeffi- 

2 

cients of the binomial expansion for a fractional power. 
Then replacing the above equation in Eq.(3.14) and 
using the binomial formula to expand ^^(d — 2) w 
[7 - (l/(d - 3)] 1/2 we obtain, for d = 3, 



—2 2 24 * 1 1 

772 ~ 777 H ^A — 7 f 



Li 



(3.17) 



Taking to 2 = a(T — T ), with a > we have from the 
above equation for to 2 = 0, the critical temperature as a 
function of the film thickness L and of the bulk transition 
temperature, T , 



T 24011 I 



5 V7 



(3.18) 



where, in agreement with the standard notation, we have 
introduced (3 = A, the Ginzburg-Landau parameter. This 
is an equation that describes the behaviour of the critical 
temperature of a superconducting film of thickness L, 
taking into account the gauge fluctuations. Of course 
when L — > 00, we recover the case of the material in bulk 
form. 

We see clearly two separated contributions in the ex- 
pression to the critical temperature given in Eq.(3.18). 
The first one due to the self-interaction of the scalar field, 
and the other coming from the interaction between the 
scalar and gauge fields. This last one would characterize 
a charged phase transition. It should be noticed that the 
self interaction contribution to the critical temperature 
depends on the inverse of the film thickness, while the 
charged contribution goes with the inverse of the square 
root of L. Taking T c = in Eq.(3.18), we obtain a posi- 
tive solution for L, 



L (0) = 



V7 



+ 



96/3 7 



Tc 



a<-2)2y/n 



(3.19) 



IV. CONCLUSIONS 

In this paper we have considered the Ginzburg-Landau 
model, confined between two parallel planes, and in the 
transverse unitarity gauge, as a model to describe a su- 
perconducting film. To generate the contributions from 
gauge fluctuations, we have used the Gaussian effective 
potential, which allows to obtain a gap equation that 
can be treated with the method of recent developments 
[13,14]. We have derived a critical equation that de- 
scribes the changes in the critical temperature due to 
confinement. Independent contributions from the self in- 
teraction of the scalar field and from the gauge field fluc- 
tuations are found. Our approach suggests a minimal 
film thickness for the existence of both charged and non- 
charged superconducting transitions. The behaviour de- 
scribed in Eqs.(3.18) and (3.19) may be contrasted with 
the linear decreasing of T c with the inverse of the film 
thickness, that has been found experimentally in materi- 
als containing transition metals, for example, in Nb [21], 
in W-Re alloys [22] and in epitaxial MgB2 films [23]; for 
some of these cases, this behaviour has been explained in 
terms of proximity, localization and Coulomb-interaction 
effects. With our formalism such a linear decreasing of 
T c with the inverse film thickness, can be directly ob- 
tained from our Eq.(3.18), simply taking e = (this could 
mean that the transitions observed in the above quoted 
references are non-charged). Also a comparison may be 
done with recent theoretical results for type II super- 
conductors [24] , where a similar behaviour of the critical 
temperature with the film thickness has been found for 
non-charged transitions. Moreover we would like to em- 
phasize, that our results do not depend on microscopic 
details of the material involved nor account for the in- 
fluence of manufacturing aspects, like the kind of sub- 
strate on which the film is deposited. In other words, 
our results emerge solely as a property that appears in 
the context of the Gaussian effective potential formal- 
ism and as a topological effect of the compactification of 
the Ginzburg-Landau model in one direction. Finnaly 
we would like to remark that our approach includes both 
charged and non-charged transitions. Our results are in 
qualitative agreement with the experimentally observed 
behaviours mentioned above, under the assumption that 
they correspond to non-charged phase transitions. 
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For L < Z/(°), the critical temperature (in absolute units) 
becomes negative, meaning that is the minimal phys- 
ically allowed film thickness, below which the supercon- 
ducting transition is suppressed. 
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